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Abstract. In a recent paper published in this journal [42 (2009) 495004] we studied a 
one-dimensional particles system where nearest particles attract with a force inversely 
proportional to a power a of their distance and coalesce upon encounter. Numerics 
yielded a distribution function h(z) for the gap between neighbouring particles, with 
h(z) ~ for small z and fi{a) > a. We can now prove analytically that in the 

strict limit of z — >• 0, j3 = a for a > 0, corresponding to the mean-field result, and we 
compute the length scale where mean-field breaks down. More generally, in that same 
limit correlations are negligible for any similar reaction model where attractive forces 
diverge with vanishing distance. The actual meaning of the measured exponent /3(a) 
remains an open question. 
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1. Introduction 

In a recent publication [T] we have studied an infinite system of particles on the 
line, located at {xi(t)}°Z_ OQ , when nearest particles attract one another with a force 
inversely proportional to a power a of the distance and coalesce upon encounter. In the 
overdamped limit, we can write 

dcCj j\. ( 1 1 \ , „ 

(i) 



(2) 



dt 

and the gaps between particles, U = Xi+i — Xi, obey 
dl A f 1 2 1 



dt a \l?_ x If lf +1 



Reacting particles systems such as this, evolving deterministically, occur often enough 
to merit further study [21 EH IH El El IS]- The case of a = 1 occurs in the study of 
crystal growth [2]. For a more detailed discussion see Ref. pQ. 

In the following we will consider the case of a > 0, which corresponds to attractive 
forces that diverge as the gap between nearest particles vanishes. The system is unstable 
and neigbouring particles tend to attract and coalesce, which leads to a reduction in the 
number of particles, i.e., to a coarsening process. 

Dimensional analysis, as well as a scaling hypothesis applied to the pertinent 
Fokker-Plank equations, yield the coarsening law [1] for the average distance i between 
nearest neighbour particles, 

lit) = [2A{l + a)t] n , n= 1/(1 + a). (3) 

This firm theoretical result is also strongly supported by unambiguous numerical 
evidence. In contrast, the distribution function h(z), for the reduced gap z = l/£ 
between particles, has proved to be more challenging, both numerically and analytically. 
In [I] we have presented numerical simulation results supporting a power-law behavior 
for small z, 

h(z) ~ z^ a) , (4) 

where (3(a) > a: (3(0) ~ 1/2, and (3(a) — > a as a — > oo (Fig. [[]). 

Analytically, we were unable to do better than mean-field, which predicts (3 = a. 
However, qualitative considerations had suggested that (3 might be larger than a, in the 
presence of sufficiently strong anti-correlations between adjacent gaps. 

Below, we show analytically that the exponent (3 and the coarsening exponent n 
are related to one another. It follows from this relationship that in the strict limit of 
z — » 0, (3 = a. The remainder of the paper is devoted to the question of why numerics 
yield a different value for (3. 

2. Relation between the exponents (3 and n 



Assuming that a > 0, Eq. (j2J) tells us that once the gap ij is small enough its own size 
dominates its rate of shrinking and ultimately this rate becomes independent of the 
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Figure 1: The exponent (3, that describes the small gap behaviour of the gap distribution 
function, h(z) ~ z" , as a function of the force exponent a, in our deterministic reaction 
model. The solid line indicates the strict z — > limit /3 — a of Eq. (II ip . After Fig. 1 of 
Ref. pQ. 



neighbouring gaps, li±\. Note that h(z) -> as 2 -> 0, so the probability that either of 
the U±i be also small can be safely neglected. Thus, in this limit, 

dk 2A 

~dt ~ ~alf ' 
whose integrations gives 



(5) 



If r is the time that it takes for a small gap of size /«(0) = 5 to close up, U{t) 
(resulting in a coalescence event), we get 



8 



a+l 



■2A[2±±,t. 

a 



(6) 
= 

(7) 



Suppose that at time t we have iV particles. Between t and t + r all gaps smaller 
than 5 would coalesce, so the fraction of coalescence events during the time r is equal 
to the fraction of intervals smaller than 5, 



dN 
~~N 



s/e 



dz h(z) 



B 



P + l 



B 



P + 1 



2A 



a 



a 



0+1 



2+1 

n + l 



(8) 
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where we have assumed the small- z behaviour h(z) = Bz^ and we have used Eq. (J7J). 
On the other hand, since £ = L/N, where L is the constant, total length of the system, 
dl = —£dN/N, or dN/N = —dlji. Differentiating Eq. ((3]) for an infinitesimal time 
dt = r, we get 

dN 2A 

(9) 



We can now equate the right hand sides of Eqs. (jS]) and ([9 

i+i 

In -+- 1 \ 

2A 



B 1 



P + 1 



a 



2A T, (10) 



which must be equal for any small r. This implies 

(3 = a and B = a. (11) 

It is worth mentioning that the exponents n and (3 can be related, in principle, 
following the same procedure as above, also when the attractive force does not diverge 
with vanishing gap (as for, e.g., an exponential force, f(l) ~ e~ l ), or when the force is 
repulsive (our model, with a < 0). However, in such cases adjacent gaps influence the 
outcome and correlations are important even in the limit of vanishing z, which precludes 
a simple derivation of S(t). 

3. Correlation effects: discussion and conclusions 

Numerically, the observed small-gap exponent, /3, seems to be larger than a, the value 
predicted from the strict limit of z — > (Fig. [T]). How small need z be for the result 
— a to hold, even approximately? 

Assuming that the small gap 5 is surrounded by gaps of typical length £, according 
to Eq. (T5]) one can neglect their influence when l/5 a ^> l/£ a , or 

j) <e, (12) 

where e is some small positive number representing our error tolerance. We can expect 
that correlations with adjacent intervals are negligible at (reduced) distances z < e l l a . 
This is consistent with the findings of Fig ([1]) that the error is larger for smaller a. 

We can make a more accurate assessment of the influence of correlations from 
neighbouring gaps. Let ^2(^1,^2) be the joint probability density for adjacent gaps of 
reduced lengths z\ and z 2 , and denote by 

i\ = r^i^ (13) 

the conditional average of z^ a , given that the adjacent interval has length z 2 = z. This 
conditional average obeys the exact expression (Eq. (30), Ref. [T]) 



2ah(z) + azh'(z) 



±) - i ) »w 



(14) 
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where the prime denotes differentiation with respect to z. Numerically, one measures 
the exponent /3 from the slope of a log- log plot of h(z) vs. z. Rearranging (I14p. we 
obtain for the local slope: 



. . d(\nh) zh(z)' 
9a(z) = — — r = -r-r-r- = 



-2a + az~ a ~ x + (z± a )' z 



(15) 



d(\nz) h(z) ' olz + z- a - (j&f a ) 

Further progress depends on At the simplest level, mean- field says that 

^2(^1,^2) = h(zi)h(z2), leading to a constant value of the conditional average, which 
can be shown to be [lj 

1\ = l + a. (16) 

Z l I z 

Another tractable possibility is to assume that adjacent gaps are perfectly anti- 
correlated: as one interval grows the adjacent gap shrinks, and their total length is 
fixed, h + h = 2£, such that h 2 (zi, z 2 ) = h(zi)6(2 — z\ — z%). This leads to 

1 X 1 (17) 



(2-zY 




Figure 2: (a) Plot of the local slope g a {z) for a = 0.001, 0.1, and 1 (inset), for the mean- 
field case (dashed lines) and for perfect correlations (solid lines), (b) The boundary z* 
within which the mean-field local slope deviates from a by less than 10 %, as obtained 
numerically from (fT5]) and f JT6|) (symbols) compared to the first-order approximation 
(II 8p (dashed line). Note the logarithmic scale. 



In Fig. |2h we plot the local slope g(z) for the two assumptions, f lT6|) and f fT7|) . 
for a = 0.001, 0.1, and 1 (inset). In all cases, it is easy to confirm analytically that 
g(0) = a, in agreement with ([IT]) . Also evident from the figure, is the fact that g(z) 
rises sharply near the origin (with infinite slope, for a < 1), which may explain how 
numerically one may observe an effective exponent larger than a. 

We can also use the expression for the local slope, ( 1T5|) . in conjunction with ( fT6|) . 
to determine the condition for a < (3 < a(l+e). Solving to first order in e we find that 



< z < z*, with 



l/a 
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(and a similar expression for perfect anti-correlations). Note the similarity to the 
criterion for neglecting adjacent intervals, derived in the previous section. In Fig. [2b we 
plot for e = 0.1, as a function of a. The conclusion is that a change in /3 as big as 
10 % occurs, for most values of a, within a very small range of z. For example, z* < 10 -3 
for a = 0.5, while the smallest gaps we could analyse using our best numerical data 
were around z ~ 0.05. Clearly, under these conditions one cannot expect to measure 
the predicted — a. 

Ultimately, however, the question remains largely open. It is possible, in principle, 
that when the true correlations are taken into account there exists a fairly wide region 
of z > z* where the local slope g(z) « > a is nearly constant. In that case, the region 
< z < z* would be analogous to a boundary layer within which the mean-field result 
prevails and beyond which a different behaviour sets in. 

Acknowledgments 

The authors thank financial support from the NSF. PP is grateful for the warm 
hospitality while visiting Potsdam during this project, and also acknowledges financial 
support from MIUR (PRIN 2007JHLPEZ). 

References 

[1] D. ben-Avraham, O. Gromenko and P. Politi, J. Phys. A: Math. Theor. 42, 495004 (2009) 
[2] P. Politi and D. ben-Avraham, Physica D 238, 156 (2009). 

[3] K. Kawasaki, M. C. Yalabik, and J. D. Gunton, Phys. Rev. B 17, 455 (1978); T. Ohta, D. Jasnow, 

and K. Kawasaki, Phys. Rev. Lett. 49, 1223 (1982). 
[4] I. Ispolatov and P. L. Krapivsky, Phys. Rev. E 53, 3154 (1996). 
[5] A. D. Rutenberg and A. J. Bray, Phys. Rev. E 50, 1900 (1994). 
[6] M. Rost and J. Krug, Physica D 88, 1 (1995). 

[7] See, for example, S. Redner, in Nonequilibrium Statistical Mechanics in One Dimension, 

V. Privman, ed., (Cambridge University Press, 1997), and references therein. 
[8] B. Derrida, Godreche, and Yekutieli, Phys. Rev. A 44, 6241 (1991). 



